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(D ■ Abstract 

00 . 

Neutrino decay has been proposed as a possible solution to the atmospheric 

T-H ■ 

^^ ■ neutrino anomaly, in the light of the recent data from the Super-Kamiokande 

^vq . experiment. We investigate this hypothesis by means of a quantitative analy 

^D . sis of the zenith angle distributions of neutrino events in Super-Kamiokande, 

o\ 

- including the latest (45 kTy) data. We find that the neutrino decay hypothesis 

^i fails to reproduce the observed distributions of muons. 
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The Super-Kamiokande (SK) experiment has confirmed, with high statistical significance, 
the anomalous fiavor composition of the atmospheric neutrino fiux. Such anomaly is found 
in all SK data samples, including (in order of increasing energy) sub-GeV e-like and /i-like 
events (SGe and SG/i) |I|, multi-GeV e-like and /x-like events (MGe and MG/x) [^, and 
upward through-going muons (UP/i) ^. In particular, all muon event samples (SG/i, MG/i, 
and UP/i) show significant distortions of the observed zenith angle distributions, as compared 
with standard expectations. The recent muon data from the MACRO 0] and Soudan-2 ||5| 
experiments, as well as from the finalized Kamiokande sample |P], are also consistent with 
the Super-Kamiokande data. 

The observed dependence of the muon deficit on both energy and direction can be beau- 
tifully explained via neutrino fiavor oscillations in the z/^ — * z/,- channel 0. Transitions into 
sterile states are also consistent with the data [§,0, as well as subdominant oscillations in 
the Ufj^ — i> z/e channel [|10|]. Determing the fiavor (s) of the oscillating partner (s) of the muon 
neutrino will represent a crucial test of such explanation(s). In the meantime, it is useful to 
challenge the oscillation hypothesis and to investigate possible alternative scenarios |TT[ . 



Neutrino decay has been recently proposed |T^ as a possible solution to the atmospheric 
neutrino anomaly. In a nutshell, the muon neutrino z/^ is assumed to have an unstable 
(decaying) component z/^, 

cos^ = (z/^|z/d) 7^ , (1) 

with mass and lifetime m^ and r^, respectively. In the parameter range of interest, the z/^ 



survival probability reads |12 



P^^ ~ sin^ ^ + cos^ ^ exp{-aL/E) , (2) 

where a = rria/Ta, and L/E is the ratio between the neutrino pathlength and energy. The 



so 



possible unstable component of z/g is experimentally constrained to be very small []T2|,|T3 
one can take Pee — 1- 

For large values of cos^ and for a ~ 0{D^/1 GeV) {D^ = 12, 800 km), the exponential 
term in Eq. (|^) can produce a detectable modulation of the muon-like event distributions 



I4 • In particular, the expected modulation seems to be roughly in agreement |jT2|] with the 
reconstructed L/ E distribution of contained SK events ||^. However, the L/ E distribution 
is not really suited to quantitative tests, since it is affected by relatively large uncertainties, 
implicit in backtracing the (unobservable) parent neutrino momentum vector from the (ob- 
served) final lepton momentum. Moreover, the L/ E distribution includes only a fraction of 
the data (the fully contained events), and mixes low-energy and high-energy events, making 
it difficult to judge how the separate data samples (SG, MG, and UP) are fitted by the decay 
solution. Therefore, we think it worthwhile to test the neutrino decay hypothesis in a more 
convincing and quantitative way, by using observable quantities (the zenith distributions of 
the observed leptons) rather than unobservable, indirect parameters such as L/ E. 

To this purpose we use, as described in [|10[ , five zenith-angle distributions of neutrino- 



induced lepton events in Super-Kamiokande, namely, SGe and SG/i (5+5 bins), MGe and 
MG/i (5+5 bins), and UP/i (10 bins), for a total of 30 data points. The data refer to the 
preliminary 45 kiloton-year sample of Super-Kamiokande, as taken from PJ^. The theoret- 
ical calculations are performed with the same technique as in [|r^, but the fiavor survival 
probability refer now to v decay [Eq. (Q)] rather than v oscillations. As in |TD[, conservative 



errors are assumed not only for the overall normalization of the expected distributions, but 
also for their shape distortions. The (dis) agreement between data and theory is quantified 
through a x^ statistic, which takes into account the strong correlations between systematics. 
Figure 1 shows the results of our x^ analysis in the plane (cos^,a), with a given in 
unit of GeV/D®. The regions at 90 and 99% C.L. are defined by x^ ~ Xmin = 4.61 and 
9.21, respectively. As quahtatively expected, the data prefer a ~ 1 and large cos^, in 
order to produce a large suppression of i/^'s coming from below. However, the absolute 
X^ is always much higher than the number of degrees of freedom, N£,p = 28 (30 data 
points — 2 free parameters). In fact, it is Xmin — §6.2 at the best-fit point [reached for 
(cos,^,a;) = (0.95,0.90)]. The very poor global fit indicates that the zenith distributions 
cannot be accounted for by neutrino decay, contrarily to the claim of [0, which was based 
on reduced and indirect data (the L/E distribution). This situation should be contrasted 



with the z/^j '^ Ur oscillation hypothesis, which gives Xmin/^DF ~ 1 both in two-flavor and 
three-flavor scenarios [0. Even using "older" Super-Kamiokande data (i.e., the published 



33 kTy sample |1|-§|), we get Xmin ~ 65 for the v decay flt, still much higher than A^df- 

Basically, neutrino decay fails to reproduce the SK zenith distributions for the following 
reason: The lower the energy, the faster the decay, the stronger the muon deflcit — a pattern 
not supported by the data. This can be better appreciated in Fig. 2, which shows the SK 
data (45 kTy) and the expectations (at the "best-flt" point) for the flve zenith angle {9) 
distributions considered in our analysis. In each bin, both the observed and the expected 
lepton rates R are normalized to the standard values in the absence of decay i?o, so that 
"no decay" corresponds to R/Ro = 1. The data are shown as dots with la error bars, while 
the decay predictions are shown as solid lines. The predictions are affected by strongly 
correlated errors (not shown), as discussed in Appendix B of [p!0| . 



In Fig. 2, the muon data show signiflcant deviations from the reference baseline R/Rq = 
1, most notably for MGfi events. The z/ decay predictions also show some (milder) deviations, 
but their agreement with the data is poor, both in normalization and in shape. Neutrino 
decay implies a muon deflcit decreasing with energy, at variance with the fact that the 
overall (6'-averaged) deflcit is about the same for both SG/i's and MG/x's (—30%). The 
shape distortions of the muon zenith distributions are also expected to be exponentially 
weaker at higher energies (i.e., for slower neutrino decay). This makes it impossible to flt at 
the same time the distorted muon distributions observed at low energy (SGyu) and at high 
energy (UP/i), and to get a strong up-down asymmetry at intermediate energy (MG/i) as 
well. Notice in Fig. 2 that, although the predicted shape of the zenith distribution appears 
to be in qualitative agreement with the data pattern for SG/i's, it is not sufficiently up-down 
asymmetric for MG/i's, and it is definitely too flat for UP/i's, where the muon suppression 
reaches the plateau P^^j ^ c| -|- s|, in disagreement with the observed cos 6'-modulation 0. 
Finally, we remark that the electron neutrinos, despite being "spectators" in the u decay 
scenario (flat SGe and MGe distributions in Fig. 2), play a role in the flt through the 



constraints on their overall rate normalization, as in the z/^ -^ v^ oscillation case [p) 



In conclusion, we have shown quantitatively that the neutrino decay hypothesis, although 
intriguing, fails to reproduce the zenith angle distributions of the Super-Kamiokande sub- 
GeV, multi-GeV, and upgoing muon data for any value of the decay parameters a and 
cos^. Even at the "best fit" point, data and expectations differ both in total rates and in 
zenith distribution shapes. Therefore, neutrino decay (at least in its simplest form [|I^]) is 
not a viable explanation of the Super-Kamiokande observations. The strong disagreement 
between data and theory was not apparent in |]12|, presumably because the experimental 



information used there was rather reduced and indirect. 

Note added. When this work was being completed, our attention was brought to the 



recent paper ||I4[, where several scenarios — alternative to neutrino oscillations — are consid- 



ered, including neutrino decay (for cos^ = 1). The authors of [|14| find that neutrino decay 



provides a very poor fit to the Super-Kamiokande data, consistently with our conclusions. 
As far as the neutrino decay hypothesis is concerned, our work has the advantage of being 
more general (cos,^ unconstrained), more refined in the statistical approach, and updated 
with latest SK data (45 kTy) m. 
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FIG. 1. Fit to the Super-Kamiokande data (45 kTy, 30 data points) in the plane of the neutrino 
decay parameters cos^ = {vp\v^ and a = nid/Td- The solid and dotted hnes are defined by 
y? — Xmin ~ 4-61 ^-iid 9.21, corresponding to 90% and 99% C.L. for two variables. The analysis 
favors a ~ 1 CeV/D© and large cos^. However, even at the best fit point there is poor agreement 
between data and theory (Xmin/-^DF ~ 86.2/28 = 3.1), indicating that v decay is not a viable 
explanation of the Super-Kamiokande observations. 



SK (45 kTy) zenith distributions at best fit (cos f = 0.95, a = 0.90 GeV/12800 l<m) 
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FIG. 2. Zenith angle distributions of Super-Kamiokande sub-GeV e-like and ^-like events (SGe 
and SG/i), multi-GeV e-like and ^-like events (MGe and MG^), and upward-going muons (UP^). 
Data: dots with itlo" statistical error bars. Theory (v decay best fit): solid curves. In each bin, 
both theoretical and experimental rates R are normalized to their standard (no decay) expectations 
Ro. The solid curves do not appear to reproduce the muon data pattern. 



